Nodal quasi-particles in stripe ordered superconductors 
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We study the properties of a quasi-one dimensional superconductor which consists of an al- 
ternating array of two inequivalent chains. This model is a simple charicature of a striped high 
temperature superconductor, and is more generally a theoretically controllable system in which the 
superconducting state emerges from a non-Fermi liquid normal state. Even in this limit, "d-wave 
like" order parameter symmetry is natural, but the superconducting state can either have a complete 
gap in the quasi-particle spectrum, or gapless "nodal" quasiparticles. We also find circumstances in 
which antiferromagnetic order (typically incommensurate) coexists with superconductivity. 



o 
o 

Oh' 



X3 

o 



> 
o 
in 
m 
o 

o 
o 

-I— > 



I 

o 
o 



X 
S3 



A key feature of the cuprate high temperature super- 
conductors is that the "normal" state is not well de- 
scribed as a Fermi liquid. Therefore, to understand the 
physics of the transition temperature, the brilliantly suc- 
cessful BCS theory, which presupposes Q| that the nor- 
mal state is a Fermi liquid, must be modified. Although 
the cuprate high temperature superconductors are lay- 
ered materials, self-organized one dimensional structures 
1^, or "stripes," have been widely observed, making plau- 
sible the idea that at intermediate scales these materi- 
als can be thought of as quasi-one dimensional. The 
only theoretically well understood example of a supercon- 
ducting system with a non Fermi-liquid (NFL) normal 
state is a quasi-one dimensional superconductor. Here 
the "normal state" is governed by the quantum criti- 
cal physics of a decoupled set of one dimensional elec- 
tron gases (IDEG); superconducting long-range order is 
triggered by inter-chain coupling, and accompanied by a 
crossover to higher dimensional physics ||] . 

Several salient features of the high temperature super- 
conductors are naturally understood from this viewpoint 
1^. The fact Q| that in under and optimally doped ma- 
terials, quasi-particles in the gap antinodal regions of the 
Brillouin Zone (BZ) exist only in the superconducting 
state, and have a quasi-particle weight {Z{T,x)) which 
vanishes as the transition point is approached, either as 
a function of the doped hole concentration, x, or of the 
temperature, T, can, in our opinion, only be understood 
if the normal state has no well defined quasiparticle ex- 
citions. The extraordinarily one dimensional dispersion 
apparent |^ in the quasi-particle spectrum in this region 
of the BZ is independent evidence of a quasi-ld origin 
for the NFL behavior. That the zero temperature super- 
fluid density, and with it T^, is roughly proportional to x 
is also simply understood in this way, as is the fact that 
the pairing scale, Aq, and the superconducting Tc appear 
to be distinct energy scales in the problem. 

We thus propose studying the high temperature super- 
conductors by adiabatic continuity from the quasi-one- 
dimensional limit. Conceptually, we imagine introducing 
an explicit symmetry breaking field of strength h into the 
physical Hamiltonian, so that for h large the problem is 



iterally quasi-one-dimensional, and can hence be solved 
using the powerful non-perturbative methods developed 
for the theory of the IDEG. Then, so long as there is no 
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FIG. 1. Schematic phase diagrams and qualitative RG 
flows of an array of alternating A and B type chains with 
Af > and A^!^ = Af = Af = 0. The figures represent cuts 
through a multidimensional parameter space. Interactions 
that couple neighboring A and B chains increase along the x 
axis, while the A to A and B to B couplings increase along 
y. The RG flows in the neighborhood of the flxed points (ex- 
cpect C4) follow from the analysis presented in the text; the 
phase boundaries and global flows are qualitative renderings. 
The unstable fixed points represent decoupled chains (a), the 
two fluid state with a 2D gapped superconductor on A and an 
anisotropic Fermi liquid on B (/?), and various critical points 
{Cj). The phases controlled by the stable flxed points are: 
coexisting superconducting and antiferromagnetic order with 
a full quasi-particle gap {13'), a 2D nodal superconductor (7), 
and a 2D fully gapped superconductor (S). 

phase transition as a function of h, the results should 
be qualitatively correct, even as /i — > 0. Moreveover, so 
long as the isotropic system has substantial local stripe 
order, many of these results should even be quantita- 
tively reasonable. Empirical evidence which suggests 
that such explicit symmetry breaking is innocuous comes 
from experiments in the strongly orthorhombic materials, 
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YBa2Cu307-5 and Y2Ba4Cu807_5, where anisotropies 
in the in-plane superfluid density as large as a factor of 
10 can be induced without, apparently, affecting the qual- 
itative physics of high temperature superconductivity. 

However, there is also experimental evidence of gap- 
less "nodal" excitations at low temperatures, deep in the 
superconducting state. In all the simplest realizations of 
a quasi-one dimensional superconductor, for example an 
array of weakly interconnected doped two-leg i — J or 
Hubbard ladders, the superconducting state has a gap 
to all spin-carrying excitations, including quasi-particles. 
The gapless and fully gapped phases of a superconduc- 
tor are necessarily separated by a quantum phase tran- 
sition! In order to substantiate the claim that the quasi- 
one-dimensional limit is adiabatically connected to the 
physics of the cuprates, it is at least necessary to de- 
termine whether, and under what circumstances, gapless 
nodal excitations exist. 

It is important to stress that the issue of the "d-wave- 
like" character of the order parameter (i.e. whether the 
expectation value of the pair creation operator changes 
sign under a 90° rotation) is distinct from the issue of the 
existence of nodal quasi-particles. Even in a weak cou- 
pling (BCS) superconductor, if the Fermi surface does 
not close around the k = (0, 0) or k — (tt, tt) points in 
the BZ, it is possible to have perfect d-wave symmetry 
without nodal quasi-particles. For the two-leg ladder of 
the above cited example, the superconducting pairing is 
known ||^ to be d-wave-like but fully gapped. 

We will study the simple model effective Hamiltonian, 
H = H* + H' , of a quasi-one dimensional system which 
consists of an array of alternating, inequivalent (A and 
B type) "chains." For zero interchain coupling, H' — 
{i.e. h — + oo), the problem is solved exactly using stan- 
dard bosonization methods; we refer to this limit as the 
decoupled fixed point, although in reality it is a high- 
dimensional manifold of fixed points, parametrized by 
the various Luttinger exponents of the IDEG's. It is 
indicated by the point a in Fig. |l|. By studying the per- 
turbative renormalization-group (RG) flows in the vicin- 
ity of the decoupled fixed point, as well as the behavior 
at various anisotropic 2D fixed points, we determine the 
qualitative phase diagram and the character of the ex- 
citations in the limit of small, but non-zero interchain 
coupling. 

The decoupled fixed point Hamiltonian H* , consists of 
a sum of terms for each decoupled chain. This problem 
can be solved exactly using methods of bosonization |^ to 
express the Hamiltonian for each chain as a sine-Gordon 
field theory for the spin and charge degrees of freedom, 
respectively. The electronic field operators on each chain 
can be expressed in terms of the charge and spin fields, 
(j)n.s and (j)n,c, and their duals, 9n,s and 9n.c^ as 

-^n,a{x) = e*'^'^^«'„,^,+ (a;) + e~''^^^*„,<,,_(x) 



*....±(X) = Un.a e'VV^[«"+'^«r±^"±-^?l (1) 

where the Klein factors obey {Afn,a,-^n' ~ 
Snn'Satr' /tto, fl is the short-distaucc cutoff, a — ±1 for up 
and down spins, kp is the Fermi wave number of chain 
n, and the bosonic fields satisfy [4>n,a{y),dxOm,a'ix)] = 
i^n,m5a,a'5{x — y). For n even, the A chains, which repre- 
sent the stripe regions where the concentration of mobile 
"doped" charges is high, are metallic but with a spin-gap 
aI'^'' > 0, i. e. a Luther-Emery liquid. Microscopically, 
one can imagine that each A chain represents the low en- 
ergy physics of a doped two- leg or three-leg t — J oi Hub- 
bard ladder ||]. The Fermi wave- number, k^j^^ is there- 
fore far from the commensurate value, 7r/2; experiment 
suggests that in a variety of cuprate superconductors, 
fc^' « 7r/4. The remaining quantities which characterize 
the A chains are a charge Luttinger parameter jl^, Ki'^\ 

[A) [A) 

and the charge and spin velocities, Vc and Vg . The B 

chains represent the more lightly doped, locally antifer- 

romagnetic strips between stripes. We will be interested 

in the case in which some doped holes have leaked into 

these strips, so they have no charge gap, but because 

they are still nearly Mott insulating |jll[, fc^"* w 7r/2, 

iiTc^-* « 1/2, and Vc^'^ <^ Vc^\ We will, however, also 

consider the cases in which Umklapp scattering opens a 
is) 

charge gap, Ac , on chains B, which, in addition, may or 
may not have a spin-gap. In the more interesting gapless 
case, spin rotation invariance implies that at low energies 
the spin Luttinger exponent, Kl — > 1. 

The interchain coupling, H' , typically [p^ generates in- 
teractions that are relevant in the RG sense at the decou- 
pled fixed point. Starting from a microscopic viewpoint, 
one would be tempted to take H' to consist of a single 
particle hopping term which couples each A chain to its 
nearest neighbor B chain. This interaction is manifestly 
irrelevant, both because of the presence of a spin-gap, 
and because of the mismatch in Fermi wave numbers, 
fc^-* 7^ k'"J^\ However, in the initial stages of renormal- 
ization, all imaginable local terms consistent with sym- 
metry are generated, both relevant and irrelevant. We 
will therefore skip this initial step, and directly study 
the perturbative /3-functions for potentially relevant in- 
terchain couplings. Here we assume that H' contains 
only the most relevant interactions, between first and sec- 
ond neighbor chains, for the range of Luttinger exponents 
discussed above, 

H'=Y, [dxl -ti3B^[1'L-i,.*2„+i,.+h.c.] 

— J_B_B^2n-l ■ 'S2n+1 — J^^B [A2„ A2„+l + h.c] 

+ ^AB[AL(*2„-l,T*2«+l,i + *2„+l,T*2„-l,i) + h.c.] 

-jAA[Ai^A2n+2+h.C.]\ (2) 
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where A„ = '^n,-\'^n,i is the singlet pair creation oper- if we set K, 



'B) 



1/2, this is true so long as K, 



{A) 



ator, S is the spin-density operator, and it is implicitly 
understood in the above expression that any piece of the 
interaction that is rapidly oscillating (with wave-number 
2kp) is to be omitted. Here, we have neglected possibly 
relevant back-scattering interactions which could poten- 
tially promote charge-density wave formation; the basis 
for this is discussed in Ref. |l2| and below. 

The only non-standard term in H' is the term propor- 
tional to J^'j\B which removes a pair from an A chain, 
rotates it by 90°, and reinserts it across the two neigh- 
boring B chains; the sign of this term (after renormaliza- 
tion, away from the decoupled, a, fixed point) determines 
whether the superconducting order is d-like or s-like. Mi- 
croscopic calculations ^ on t — J and Hubbard ladders 
lead us to believe that most likely all the pair-tunnelling 
terms Jx > 0. The d-like pairing tendency observed in 
these calculations implies JT^^ > 0, while the positivity 
of the remaining JJx imply an unfrustrated supercon- 
ducting state; e.g. Jab < 0, would imply "7r-junctions" 
between neighboring chains. 

We begin by considering the regime in which all of the 
coupling constants in H' are small. In this limit, the 
system is at the decoupled fixed point a. It is easy to de- 
termine the (perturbative) role of the various processes 
in H' . In cases in which a gap prohibits the operation of 
one of these interactions in lowest order, that interaction 
is manifestly irrelevant, i. e. its dimension is infinite. 
Otherwise, to leading order in powers of the interchain 
couplings, the perturbative /3- functions have the form 



c?lna 



= i2-D,)g 



(3) 



Dg is the scaling dimension of the perturbation with cou- 
pling constant g. For Dg < 2, the operator is perturba- 
tively relevant, and otherwise it is irrelevant. The di- 
mension of the various operators are listed in Table |. 
Forward scattering interactions between the charge cur- 
rents and densities on neighboring chains are precisely 
marginal. 

We now return to the issue of the CDW couplings be- 
tween chains; if relevant, these would lead to an ordered, 
insulating state. In Ref. 
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it was shown that forward 
scattering interactions between chains, whether direct or 
induced by dynamical fluctuations of the stripe geometry, 
strongly affect the scaling dimension of these operators, 
tending to make them less relevant. In particular, there is 
a finite regime of parameters, especially when the stripe 
fluctuations are significant, where the CDW couplings 
are irrelevant, and so can be neglected. 

Since in all the cases considered here, the decoupled 
fixed point is perturbatively unstable in some way, our 
next task is to determine where the RG flows go. 

The two fluid fixed point: It should be clear from the 
table that, under most conditions, all AB couplings are 
perturbatively irrelevant at point a of Fig. 0. Specifically, 



< 1. 

(The phase diagrams in Figs, la, Ic, and Id assume this 
condition is satisfied.) Even if K^"^^ > 1 (as in Fig. lb), 
the AB couplings are typically more weakly relevant than 
the couplings between like chains. We are therefore led 
to consider the RG flows in the limit in which all such 
couplings are set equal to zero, so that we have two in- 
terpenetrating, and decoupled, but genuinely two dimen- 
sional systems. This limit is represented by the left-hand 
edge of the phase diagrams in Fig. ^ 

In general, there are many possible higher dimensional 
fixed points to which the system could flow. Because 
the system can lower its kinetic energy by allowing pairs 
to move between chains Jaa is typically relevant- 
technically so long as Kc > 1/2. We will consider ex- 
clusively the state in which the A subsystem has true 
superconducting long range order and a full spin-gap. 
Because tsB always has lower scaling dimension than 
Jbb so long as Af — the most likely situation is that 
system B forms a highly anisotropic Fermi liquid. The 
corresponding fixed point is labelled (3 in Figs, la-lc. 
However, either if > 0, or if the residual interactions 
are sufhcently strong to drive an instability of the nearly 
nested Fermi surface, it is possible that system B will or- 
der antiferromagnetically. The fixed point with coexist- 
ing superconductivity in A and antiferromagnetism in B 
is indicated by /3' in Fig. Id. (Quantum fluctuations can 
also lead to a magnetically disordered insulating state on 
system B, which is adiabatically connected to the state 

( B) (B) 

in which Ac and As are non-zero.) 

The operators which couple A and B have entirely dif- 
ferent scaling dimensions at the two-fluid fixed points, /? 
and /?', than at a. The two-fluid fixed point (3 is un- 
stable, due to the ordinary proximity effect. Specifically, 
in establishing the relevance of the operators Jab and 
J^g , the pair creation operator on A can be replaced by 
their constant expectation value, while the pairing fields 
on B operate on a Fermi liquid state. For Jab > this 
superconducting state is d-wave-like; the quasi-particle 
spectrum is: 
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TABLE I. Scaling dimensions of the interchain couplings 
at the decoupled fixed point with A^ > and A^ = 0. 
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E^(k)= [vf{K - kf) - 2tBB cos(fcj,L)]2 
+ <Aa>^ [Jab - '^J'ab cos{kyL)f 



(4) 



where x and y refer, respectively to the directions parallel 
and perpendicular to the chains, vp, tsB, Jab and J'ab 
are now to be interpreted as renormalized parameters at 
/9, and L is the spacing between chains. There remains 
the quantitative issue, which depends on the microscopic 
details, of the relative magnitudes of these couplings. If 
\Jab\ < %J'ab\^ the quasi-particle spectrum is gapless. 
This is a stable phase of matter, the nodal superconduc- 
tor, described by the fixed point 7 in Fig. 1. We have as- 
sumed that this inequality is satisfied in Figs, la and lb, 
where the RG flows run directly from ji to 7. However, if 
\Jab\ > '^\Jab\' flows run from (3 to the fully gapped 
superconducting state, signified by the stable fixed point 
d in the figure, as in Fig. Ic. 

By contrast, at the fixed point /3', there are no finite 
dimensional operators which couple the two fluids since 
the antiferromagnet has a charge gap and the supercon- 
ductor has a spin-gap. This fixed point describes a stable 
phase with two coexisting order parameters, and a com- 
plete gap in the quasi-particle spectrum. Evidence of 
a phase with coexisting magnetic and superconducting 
order, presumably rendered glassy by quenched disorder, 
has been presented 1 13 ij] ; it is a prediction of the present 
study that the ordered state should have a fully gapped 
quasi-particle spectrum. 

Fully coupled fixed points: In addition to the stable 
fixed points 7 and S, described above, there are a number 
of other unstable fixed points whose existence is dictated 
by the topology of the phase diagram (assuming that the 
transitions are second order). 

The phase transition between the nodal and the 
gapped superconductor is governed by the fixed point 
Ci in the figure. The universal properties of this transi- 
tion can be studied [jisl in the weakly interacting limit. 
This transition, as approached from the nodal phase, is 
triggered by deforming the gap function (or band struc- 
ture) such that two nodal points approach each other, 
and at the critical point, coalesce. As a consequence, at 
the single nodal point, the quasi-particle velocity in one 
direction vanishes. Beyond the critical point, a gap opens 
in the spectrum. By naive power counting, four fermion 
interactions are irrelevant at this fixed point, so this is 
all there is to it! Even when Jab is perturbatively irrel- 
evant at a, strong enough pair-tunnelling will certainly 
produce superconductivity in B by the proximity effect. 
Moreover, the resulting state will be the fully gapped, 
ordered superconductor. Thus, there must exist a highly 
unstable critical point, C2, similar to the smectic metal 
to superconductor critical point of Rcf. . 

The critical point C3 in Fig. Id is similar to one 
mentioned by Vojta et. al. ]l5[ |, for a transition, within 
a nodal superconducting phase, to a state with broken 
translational symmetry with an ordering vector which at 



C3 spans the nodal points. Conversely, the robustness 
of the nodal superconducting phase embodied in Figs, 
la, lb, and Id, mirrors the asymptotic decoupling of the 
nodal quasi-particles from fluctuations associated with 
such a transition when the ordering vector does not span 
the nodes. The (as yet not analyzed) multi-critical point 
C4 is required by a minimal consistent construction of 
the RG flows. 

Finally, it is clear that at the fixed point 7 the nodal 
quasi-particles are well defined elementary excitations 
of the system at arbitrarily low energy and long wave- 
length. This is consistent with a widely held belief that 
such excitations must exist, even though the high temper- 
ature superconductors are far from the BCS limit. How- 
ever, what is also clear from the circuitous flows that 
lead to this fixed point in Figs, la and lb is that these 
quasi-particles can, under appropriate circumstances, be 
much less robust than the superconducting state, itself. 
Thus, even well below Tc, where superconducting order 
is well established, the nodal quasi-particles can still be 
ill-defined objects, and only become sharp at very low 
temperatures. It is possible that this observation recon- 
ciles the strong evidence jl6| from thermal conductivity 
of well defined nodal quasi-particles at T <C Tc with the 
evidence from photo-emission |l7|] of their non-existence 
down to temperatures of order Tc/2. 
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